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Abstract 

An asymptotic theory is developed for a moving drop driven by a wettability gradient. We distinguish the 
mesoscale where an exact solution is known for the properly simplified problem. This solution is matched 
at both - the advancing and the receding side - to respective solutions of the problem on the microscale. 
On the microscale the velocity of movement is used as the small parameter of an asymptotic expansion. 
Matching gives the droplet shape, velocity of movement as a function of the imposed wettability gradient 
and droplet volume. 
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I. INTRODUCTION 



The description of the movement of a three-phase contact line is an essentially unsolved hydro- 
dynamical problem that continues to attract much interest, for instance, when studying spreading 
drops, and liquid sheets or ridges moving down an inclined plate. The understanding of 'simple' 
contact line movement is also paramount for a deeper insight in related problems as the dynamical 
wetting transition and transversal instabilities of moving contact lines. 

It is well known that the divergent shear stress at the contact line forbids a solution in the 
framework of purely classical hydrodynamics, i.e. assuming a no-slip boundary condition at the 
solid-liquid interface. Although this was first pointed out by Huh and Scriven 1 based on Moffatt's 2 
solution for flow in the edge which does not satisfy the normal stress boundary condition on a free 
interface, and was never proven rigorously, the divergence can be understood as a consequence 
of incompatibility of multivaluedness of the velocity at the contact point in the classical hydrody- 
namic formulation. 

The boundary condition has to be relaxed to permit movement of the contact line. This can be 
done by introducing a very thin precursor film on the 'dry' substrate 3 , or by allowing for slip at 
the solid-liquid interface everywhere 1 or only near the contact line 4 ^, or introducing an effective 
molecular interaction between the substrate and liquid into the hydrodynamic model 6 -. For a 
discussion of the slip condition see also the review by Dussan 8 . Other approaches include phase 
changes at the contact line 9 or introduce the vapour-liquid or fluid-solid interface, or both, as 
separate phases with properties that differ from the bulk fluid 10 . 

Most of the work on moving liquid sheets and ridges prescribes a precursor film or slip at the 
substrate. Divergence problems at the contact line are avoided, but at the expense of introducing 
ad hoc parameters into the theory. These, namely the slip length or the precursor film thick- 
ness, influence the profile of ridges and fronts and hence also the characteristics of the transverse 
instabilit y 3 ' 11 , 12 ! 13 . 

The most realistic option is the explicit introduction of molecular interactions into the hydro- 
dynamic formalism. This is accomplished by means of an additional pressure term, the disjoining 
pressure 14 . Depending on the particular problem treated, this disjoining pressure may incorpo- 
rate long-range van der Waals and/or various types of short-range interaction term o 6 ' 7 ' 15 . Recently 
Pismen 16 derived a film thickness equation with a disjoining pressure term by combining the long 
wave approximation for thin films 17 with a nonlocal diffuse interface description for the liquid-gas 
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interface that incorporates van der Waals interactions. 

These interactions are essential for the process of dewetting, and studies of dewetting 
of a thin liquid film on a substrate are generally based on models involving a disjoining 
pressure 18 i 19 i 20 i 21 i 22 i 23 . Only a few studies of instabilities of liquid fronts have adopted a similar 
approac h 24 ' 25 ' 26 , despite the fact that such an approach predicts all the ad hoc parameters of the 
slip or precursor models (i.e., the static and dynamic contact angle, drop velocity, and the drop 
and precursor film thickness) connected with the wetting properties of the liquid in terms of the 
parameters characterizing the disjoining pressure. 

Recently, Eggers presented asymptotic solutions for the profile of advancing 27 and receding 28 
driven contact lines (see also Ref.l^l). The respective solutions match inner solutions near the 
contact line where a slip model is used and outer solutions based on an analytic solution in terms 



of Airy functions discussed in Refs. 



However, the advancing and receding case are 



studied for a plate pushed into and pulled out of a liquid bath, respectively. It is not possible to 
directly couple the two asymptotic solutions to describe the motion of a driven moving droplet or 
ridge. 

In the present work we tackle the problem of an asymptotic description of a gradient-driven 
moving droplet that encompasses both an advancing and a receding contact line. This implies 
that the description of the two contact lines and the respective matching procedures depend on 
each other. Thereby we explicitly introduce the molecular interactions into the hydrodynamic 
formalism by using a chemical potential or disjoining pressure describing a situation of partial 
wetting. This corresponds to a precursor film model where the precursor film thickness is deter- 
mined through the disjoining pressure. 

We distinguish among three regions: 

• Microscopic (molecular scale) region: the dominant balance is between disjoining potential 
and surface tension. 

• Mesoscopic region: the dominant balance is between viscous dissipation and surface ten- 
sion. 

• Macroscopic region: the dominant balance is between surface tension and external forces 

Examples for driving forces are gravity for droplets or fronts on inclined plates, Marangoni forces 
occurring if temperature gradients along the substrate exist or wettability gradients along the sub- 
strate. Both, gravity and Marangoni forces act in the lubrication limit as bulk forces, i.e. the force 



is fed into the system in a top-down manner. This implies that the macroscopic region has to be 
included in the description. However, the third mentioned way to drive the system is based on 
a force resulting from a wettability gradient that is fed into the system in a bottom-up manner, 
i.e. on the microscale. The simplest description of such a system is undertaken here by matching 
solutions obtained in the mesoscopic and microscopic region. If the droplets are small enough 
(smaller than the capillary length) the macroscopic scale can be ignored. 

There are different physical situations where a gradient in wettability occurs that can be mapped 
onto the presently studied model, (i) A droplet can 'sit' on a step in wettability 34 allowing for 
an intermittent range of stationary movement until the complete drop sits on the more wettable 
substrate, (ii) A droplet can move along a smooth wettability gradien t 35 ' 36 ' 37 ' 38 , (iii) In a situation 
involving an adsorption reaction at the substrate underneath the droplet a droplet can produce the 
wettability gradient that drives its mo vemen t 39 ' 40 ' 4 1 ' 42 . In this way it carries the gradient along 
with its movement. The latter case is also related to droplet motion caused by a surface phase 
transitions^ 3 -. 

In the following we study all these situations in a model that uses a chemical potential with 
different constants at the advancing and the receding contact line, respectively. For situations (ii) 
and (iii) this corresponds to the assumption that the wettability gradient is small as compared to 
the size of the contact zone but sizable as compared to the overall droplet size. 

In the next section the basic equations for the lubrication description of moving droplets are 
introduced. The exact solution in the mesoscopic region and its asymptotics are described in 
Section [Oil The microscopic solution and its asymptotic matching are discussed in Section Hvl 
Finally, a comparison of asymptotic and numerical results is given together with our conclusions 
in section fVl 

II. BASIC EQUATIONS 

Our starting point is the thin film evolution equation in lubrication approximation 



Here 7 is the surface tension of the liquid; e is a scale ratio used as a small parameter of the lubri- 
cation expansion (which will further be identified with the local equilibrium contact angle). We 
shall use the simplest mobility function k(h) = ?7 — 1 /i 3 /3, obtained under assumption of constant 




(1) 
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dynamic viscosity 77 with no slip at the substrate. The chemical potential /j s accounts for wetting 
properties. Note, that it corresponds to the negative of a disjoining pressure II as used, for instance, 
in Ref.l^ For specific computations, we shall use the form 16 



%{h) = ^ 



(2) 



where Q s is a characteristic excess fluid-substrate interaction energy, which is proportional to 
the Hamaker constant^. If Q s > 0, this form corresponds to a negative long-range and positive 
short-range part of the spreading coefficient, thereby combining a destabilizing long-range and 
a stabilizing short-range van der Waals interaction. The contact angle is finite, and bulk fluid 
coexists at II = 0, i.e. in a flat layer of macroscopic thickness in the absence of external forces, 
with an ultrathin precursor of thickness h m . 

The variables in Eq. © are still dimensional but scaled to conform with the lubrication approx- 
imation. They are related to the physical variables (marked by a hat) as follows: 

h = h, x = x/e, t = t/e 2 . (3) 

In consequence the scaled contact angle 9 is related to the physical one by 9 = 9/e; the scaled 
droplet volume V = J hdx is related to the physical one V = J hdx by V = eV. Without any 
gradient parallel to the substrate, this model describes droplets with a finite equilibrium contact 
angle sitting on an ultrathin precursor film. 

However, here we are interested in moving droplets driven by wettability gradients along the 
substrate. In the chemical potential chosen here [Eq. ©] a wettability increase can be modelled by 
a decrease of Q s or by an increase of h m . We chose here the former possibility. Note, however, that 
in a real physical system both parameters are affected. The analysis then involves more algebra 
but is also straightforward. 

We shall consider stationary motion of a 2D droplet with the velocity U . Replacing in Eq. CO) 
d t h by —Uhxje and integrating once yields, after dropping the bars, 

S(h — h m ) d 



h 3 dx 



[h"(x)-fi s (h)} (4) 



x 3Ca r]U p, a (h) 

with b = — —, Ca = — , ^s\h) 



e A 7 



where 5 is the appropriately rescaled capillary number Ca. For 5 1, this equation is solved 
separately in the microscopic and mesoscopic regions, and solutions are matched considering a 



respective subdominant term as a perturbation. Since Eq. © does not contain the coordinate 
explicitly, the order can be further reduced (for a monotonic section) by replacing the variable 

y(h) = [h'{x)f: 

± 5{h Jyv ] = \y"( h ) - vi h m) = o. (5) 

In the next sections, solutions are determined in the mesoscopic and microscopic region, respec- 
tively. For comparison, the stationary moving droplets described by Eq. © will also be computed 
numerically using continuation technique a 45 ' 46 employing the software AUTC^T 4 ^. 

III. EXACT MESOSCOPIC SOLUTION 
A. General solution 

At large distances (h ^> 1) a simplified "mesoscopic" equation can be obtained by discarding 
the disjoining potential term in Eq. © and neglecting also h m <C h in the viscous term: 

Sh~ 2 = d xxx h. (6) 

Rescaling the height h = 5 1 / 3 ^ reduces Eq. © to a parameterless form 

C 2 = d xxx (. (7) 

This equation is invariant to simultaneous rescaling of ( and x. We chose 5 > 0, however, results 

for 5 < can be obtained by the transformation x — > —x. 

Equation © has an exact solution expressed in a parametric form through Airy functions 32 : 

K 2 X ^K 
= 2 2 f v X ^ = [ Ai (s)Bi(s ) - Ai(s )Bi(s)] 

7T Z U Z {S) U{S) 

with u{s) = AiO)Bi'(so) - Ai'(s )Bi(s). (8) 

An indefinite factor K appears here due to scale invariance of Eq. ©. It corresponds to the 
height of the droplet expressed in units of h m , i.e. it has to be large. 

The parametric solution ®, generally, defines a discontinuous function C( x )> which is phys- 
ically relevant only within certain intervals. For s < s\ where w —1.01879 is the largest 
zero of Ai'(s), physically irrelevant solutions arise with £ — > oo for x — ► ±oo and a minimum in 
between. For < so solutions exist with h — > at s — >• oo, which correspond to a sharp receding 
contact line at 

x * = 2 1 / 3 fT Ai(s )/Ai'(s ) < 0. (9) 
6 




x/K 

FIG. 1 : Droplet profiles on the mesoscale given by Eq. ® for different > so as given in the legend. The y- 
axis represents the height as Q{x)/K, i.e. corresponding to h/KS 1 ^ 3 in the scaling used on the microscale. 
For comparison, we also give the parabolic shape of a static droplet on a homogeneous substrate. The 
droplet has the same volume J Q{x)dx as the one for sq = 2. 

For < so < the height £ increases monotonically with x. These solutions are used as a 
model for a receding contact line in Ref.|3 If> however, s > 0, the profile ((x) has a maximum 
at s = s corresponding to x = 0; the solutions pass through a minimum at s min < s , i.e. 
Xmin > 0, before diverging as £ ~ x 2 , x — >• oo at s = s*, where s*(so) < is the largest zero 
of u(s). As s increases, the minimum comes very close to the x axis and the curvature at the 
minimum becomes very large. Examples of solutions for different < s are shown in Fig[T] 

B. Physically relevant interval 

We focus here on the case of moderately large so taking the profile between £ = for x = x* 
and the minimum of ((x) at x = x min as the outer solution for a moving droplet driven by a force 
fed in on the microscale. The two parameters so and K, as well as the droplet velocity that is 
absorbed into the scaling, should be obtained by matching the two inner (microscopic) solutions 
at advancing and receding sides, as well as fixing the droplet volume. 

For moderately large s one finds that the location of the minimum closely approaches s min = 
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FIG. 2: The residual height at the minimum of the droplet profiles as a function of s - The ordinate 
represents the minimal height as Cmin/K, i.e. corresponding to h m i n / K5 1 ^ in the scaling used on the 
microscale. The solid line gives the result using the approximative s m i n = —1.01879 valid for moderately 
large sq [Eq. <15b 1. The dashed line gives the result using the full Eq. (J8J). 



— 1.01879, which is the largest zero of Ai'(s). This follows from asymptotic relations applicable 
at moderately large s that will be further discussed in Section lTlICI The resulting residual profile 
height at the minimum is plotted in Fig. [2j For a physical precursor film thickness h m of the order 
of 1 nm for a millimetric drop K = 10 6 and one needs s ~ 4.5. For a droplet of one micron 
height K = 10 3 and s « 2.5. 

Using the asymptotics of Airy functions at s — > oo, 



M(s) 



-2 3/2 _ 1/4 

e 3 s 1 



Bifsl 



2 s 3/2 _ 1/4 

es S 1 



7= +0 (* 
In v 



the solution can be expanded near the zero of £ as 

2 1 / 3 A- 



x(s) 



Ai'(ao) 



Ai(s ) 



_i 3/2 

e 3 A 
27rAi'(s ) 



l + 0(s 




^ e " |s3/2 [l + 0( S - 3 / 2 )] + 0(e^^ /2 ), 



7T [Ai'(So) 

Explicit asymptotics is obtained by solving Eq. dTTT> with respect to s: 



3 x x 

_ 7 m t — 

4 L 



* \ 2 /3 r 



1 + In 



(10) 



(ID 
(12) 



(13) 



where L 1 = 2 2 / 3 irAi'(s ) 2 /K. This yields, up to corrections of higher order in In [(a; — x*)/L], 



C(x) x (x-x*) -3 In 



■* \ 1/3 



-3 In 



*\ 1/3 



(14) 



The length L is very large for moderately large sq. 



C. Limit of weak driving 



Although Eq. © does not contain the rescaled capillary number 5, we expect the applicable 
outer solution to become symmetric, approaching a parabolic profile, in the limit 5^0, which 
corresponds to a vanishing wettability gradient. As illustrated in Fig.[T] and confirmed by the 
following asymptotic analysis, the outer solution becomes almost symmetric at large values of so, 
which, as we shall further see, correspond to small values of S. 

The limit s — ► oo can be obtained with the help of the asymptotics (TTOb of Airy functions, 
which is practically applicable already at moderately large values so > 2. The resulting asymptotic 
profile height at the minimum on the advancing edge is 



( min x K{nAi(s min )Bi'(s ) 



1. 10937/0?, 



-1/2 



4 3/2 r 



e 3 a 



l + 0(*o 8/a ) 



(15) 
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The minimum is located, up to an exponentially small correction proportional to e~3 s o , at the 
largest zero of Ai'(s), i.e. s min = —1.01879 (see Fig. 01. 

The asymptotic expression for the second derivative c = ("(x) at the minimum is 



C St 



2 1/3 K -1„2 



7T 



Ai'(s T 



,<> 'rm?i ^ i (•''mm) ^i (sq 



i + o(s - 3/2 ; 



(16) 



withe w 1.15697. 

The corresponding asymptotic value of the coordinate x is 



i + K 3/2 + °( s o a ) 



(17) 



In the leading order, this coincides by the absolute value with the asymptotics of x* given by 
Eq. ©: 



x 



-2^Ks 1/2 



i -K 3/2 + °(v 



(18) 



This points out to the symmetry that should be attained in the limit of zero velocity. The full 
profile away from the location of the minimum should be computed by assuming both s and s to 



be large. This yields, in the leading order, 

C x Kis/so) 1 / 2 sech 2 [|(so 3/2 - s^ 2 )} , x x 2 1 / 3 Ks^ 1/2 tanh [f (s 3/2 - s^ 2 )} . (19) 

At the receding edge, the asymptotics of these expressions at s ^> s coincides with the asymp- 
totics of Eqs. (fTTT) . (fT2l at large s . 

The height is of the same order of magnitude as the macroscopic length scale K only when s 
is close to so- Setting s = sq, one can see by combining the above expressions that the function 
C(x) indeed approaches in this limit the parabolic profile Q/K — 1 — b 2 x 2 with b = 2~ l l 3 s]J 2 / K 
everywhere except the immediate vicinity of both contact lines; the corrections are of 0(sq 2 ). 
Thus, the scaled droplet volume is computed as 

v = 2^/3 r (1 _ bV)dx = **£L = m^i. m 

Jo 36 Ss^ 2 

IV. MICROSCOPIC SOLUTION AND MATCHING 
A. Expansion in 8 

In the microscopic region, the thickness changes from h = h m to a "mesoscopic" value far 
exceeding h m but small compared to the drop size and capillary length (that is here infinite). 
Solving Eq. © with 5 = defines the static contact angle in the limit h — » oo, while for 5 ^ 
an apparent dynamic contact angle is obtained in this limit. The appropriate length scale in this 
region is h m ; the respective dimensionless form of Eqs. @ and © is 

6^ = 4z[^)-^(h)l (21) 



h 3 dx 
h-1 



±5 h -^= l -y"{h)-pJ s {h), (22) 



with 

P 2 (, 1 



= > ^7h7^ (23) 



'in 




To model different wettability at the advancing and receding contact line one assumes different 
constants (3 = f3 a dv and f3 = (3 rec , respectively. A higher wettability at the advancing side is 
assured by f3 adv < (3 rec . 

For the receding meniscus, the positive sign should be chosen in Eq. (|22b . and the boundary 
conditions are h = 1, h'(x) = at x — » — oo, and h"(x) — > at x — >• oo, or = at /i = 1 and 
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y'(h) — > at h — ► oo. The latter condition, suggested by Eggers 28 , should fit the curvature of the 
mesoscopic solution, which, according to Eq. (TBI) , approaches — oo in the limit x — > x*. 

For the advancing meniscus, the negative sign should be chosen. The boundary condition 
h'(x) = should be set at x — > oo, and the condition h"(x) — ► at a; — ► — oo. 

The solution of Eq. (l2"2l) is sought for as an expansion in 5: y(h) = yo(h) + 5yi(h) + . . .. The 
zero-order equation, 

^(h) - rf s {h) = 0, (24) 

is readily integrated to obtain 

yo(h) = \ f? (~ + \h + 2h* + h*) , (25) 

The equilibrium contact angle #0 is obtained from the zero-order equation (l25l) in the limit h — > 00: 

9 = h'(oo) = ^0(00) = y/3/6/3. (26) 

The formal small parameter e can now be identified with the small physical equilibrium contact 
angle, say, 9q £C and expressed through physical parameters by requiring 0Q ec /e = 6l ec = 1. This 
yields 

The latter estimate follows from the estimate for surface tension 7 oc Qi/d 2 , where Qi is a charac- 
teristic interaction energy of fluid molecules and d < h m is the nominal molecular diameter. The 
contact angle is indeed small when Q r s ec / Qi (the dimensionless Hamaker constant at the advancing 
contact line) is small. The numerical value of (3 rec is specific to the particular expression for the 
disjoining potential ©, but the general procedure would be the same for any potential of a similar 
shape. Note, that now only (3 a d v < y^5/3 determines the driving wettability gradient. 

Further derivation is carried out separately for receding and advancing menisci, in view of 
different boundary conditions for the two cases. 



B. Receding meniscus 

The first-order equation derived from Eq. d2"2"b is 

h-1 1 
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Using here Eq. d25t and integrating from h = 1 to oo yields the value of y[(l) necessary to satisfy 
the asymptotic boundary condition of vanishing curvature at h — ► oo for the receding meniscus: 



2 4 

h ( - + - h + 2h 2 + h 3 
1 3 3 



-1/2 



dh 



« -1.3383 J~- = -1.3383. (29) 

V 3 (J 

The latter value corresponds to the scaling d27t . 

A non-zero value of y[(l) appears to change qualitatively the character of decay to the equilib- 
rium precursor thickness at very small deviations h — 1 < 0(5). At these distances, the expansion, 
in fact, breaks down, but the solution can be readily found by linearizing Eq. (I2T1) near h = 1. 
The linear equation is solved by a combination of exponents e Xx where A is a positive root of 
A 3 — 3/9 2 A — 5 = 0. While for 5 = the layer thickness decays at x — > — oo to unity as e^ 35 , 
for 5 7^ an additional small root A = 5/5~ 2 /3 + 0(5 2 ) appears. This root is positive, indicating a 
very slow decay to the equilibrium precursor thickness (and, possibly, breakdown of quasistation- 
ary approximation) at 5 — > 0. 

Since y[(h) ~ h^ 1 at h — > oo, |/i(/i) diverges logarithmically in the outer limit. The asymptotic 
expression is obtained by integrating Eq. (l28l) with the boundary condition (l29l) : 

yi x - 2 In — , oi« 0.444. (30) 
ai 

The respective expansion for the slope useful for further matching to an outer solution is 

h'(x)-l-5ln— + 0(5 2 ). (31) 
ai 

The expansion can be routinely continued to higher orders with the help of a symbolic computation 
program. 

To match the mesoscopic and microscopic solutions at the receding side, we compare the outer 
limit of the receding microscopic solution (h — > oo) with the inner limit of the mesoscopic solution 
(( — > 0). This translates to comparing [C'^)] 3 given by Eq. (fT4l with that given by Eq. (OTT) . After 
rescaling Eq. (fT4l and shifting the location of the contact line to zero and rearranging Eq. (l3ll . 
this gives 

[h'(x)f = -35 In - = -35 In -. (32) 

Li d\ 

The matching requirement yields the dependence of L and, hence of s on 5, expressed in an 
implicit form 

1 



L = Ai (s ) w — — exp 
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4 3/2 
Sn 

3 



a\ exp 



35 



(33) 



The approximate expression is valid for s 3> 1 (practically, for s > 2). This expression connects 
so with the dimensionless velocity 5 as shown for different K in Fig.|3l 



to 




FIG. 3: The dependence of the reseated capillary number 5 on the parameter sq for different droplet sizes 
parametrized by K (given in the legend) as given by Eq. d33l . 

Combining this result with the dependence of s on droplet size discussed at Fig.|2] gives an 
estimation of the velocity 5. For a physical precursor film thickness h m of the order of 1 nm for 
a millimetric drop K = 10 6 , so ~ 4.5 and in consequence 5 ~ 0.01 For a droplet of one micron 
height 5 » 0.025. 

Note that this is still only an order-of-magnitude estimate, because K itself depends in a subtle 
way on the velocity. This can be seen in Fig.[4]where the dependency of the droplet volume 

V = 5 1/3 K 2 f mm ((s)x'(s)ds (34) 

J oo 

on so is plotted. To compare droplets of identical volume for different driving forces, one has to 
determine K using the matching at the advancing edge. 

C. Advancing meniscus 

For an advancing contact line, the mesoscopic solution has no logarithmic asymptotics, and 
for matching one can use the zero-order microscopic solution, matching its limit at h — > oo, 
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FIG. 4: Dependence of the droplet volume on the parameter sq. The y-axis represents the scaled volume 
py#2£i/3 with v 

given by Eq. V7#o represents a physical volume that may be used as a fixed control 
parameter. The numerical result using the minimum calculated with the full Eq. (JSJi as integration boundary 
can not be distinguished from the solid line. 

h"(x) — > to the mesoscopic solution at the inflection point C"{ x ) = 0. 

This translates to comparing S 1 ^ 3 d x C(s i ) = 5 1 ^ 3 C'(si)/x'(si) at the inflection point s = given 



to h'(x — > oo) = defined by Eq. (l26l) with /3 replaced by /3 a( fo. As result of the matching, one 
finds 



i.e. 5/(0Q dv ) 3 can be calculated as a function of the parameter s as presented in Fig.|3 

This procedure effectively cuts off the highly curved segment of the mesoscopic solution near 
the minimum. Take note that Q is still much larger than the 0(1) microscopic scale, and one can 
expect corrections due to the disjoining potential to become significant only well below this value. 
However, as we will illustrate in the Conclusion, the first order matching is already sufficient to 
completely describe the droplets driven by a wettability gradient. 



by 



Sl [Ai'(s )Bi(s t ) - Ai( Sj )Bi'(s )] 2 = [Ai'(s )Bi'( Sl ) - Ai'( Si )Bi'(so)] 2 



(35) 




(36) 
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FIG. 5: Dependence of the scaled droplet velocity 5/(0Q dv ) 3 on the parameter Sq, as obtained from the 
inflection point matching at the advancing contact zone [Eq. d36l l. 

V. CONCLUSION 

We have developed an asymptotic theory for a moving drop driven by a wettability gradient. 
Wide separation between the meso- and microscale allows us to use respective analytical and 
expanded solutions on the different scales. 

Matching of the mesoscale and microscale solutions at the advancing and the receding contact 
region allows to obtain the droplet shape and the velocity of movement as functions of the imposed 
wettability gradient and droplet volume. In this way, the two matching procedures together with 
a translation between the different scalings gives a complete characterization of the droplet mo- 
tion for a given physical volume V/9 and the wettability gradient characterized by the physical 
receding 6 r ec = 6 6 T ec = 9 and advancing 9 ldv = 9o0Q dv < equilibrium angle. 

Fixing the overall long-wave scaling by fixing 9 , the three relations between 0Q dv , V, 5, K 
and s obtained in the course of the present work [Eqs. (Till). (1331) and ( BBT ) illustrated in Figs. [21 
|U and El respectively] allow us to determine the unknown 5, K and s for each given pair of 9q 1v 
and V. In Fig. |6] results of the asymptotic matching are given for the velocity 5, and the "shape 
parameter" s in dependence of the advancing equilibrium contact angle ldv for a selection of 
volumes V. As expected, the velocity goes towards zero as the driving wettability difference 
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0.03 




FIG. 6: Results of the numerical matching proceedure described in the main text using at the advancing 
contact zone the inflection point matching. Shown are (a) the droplet velocity 5, and (b) the parameter so 
describing the mesoscopic shape in their dependency on the imposed equilibrium advancing contact angle 
Gadv < &rec = 1 for different given droplet volumes V as given in the legend. In (a) results are shown for 

2 < s < 20. 



grec _ gadv van [ s \Yes, i.e. 6 r Q ec — ► 1. The shape parameter sq diverges for 9q £C — > 1 as discussed in 
section UlI CI At a fixed driving 0Q dv , the droplet becomes more asymmetrical (s decreases) and 
faster with decreasing volume. The velocity changes with volume are more pronounced for larger 
driving (i.e. smaller 0Q dv ). 

Albeit the matching is based on an expansion in 5, the numerical calulations leading to Fig.|6] 
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are not practical for very small 5 (i.e. large s ) because already for s = 20 the calculation involves 
small numbers of the order of e~ 120 that are difficult to handle. Under these conditions, the shape, 
however, remains almost static, and the integral relations of Ref. |48| can be used. The relations 
obtained by multiplying Eq. @ by h—h m and integrating over the entire x axis yield the expression 
for the dimensionless velocity in the form of a ratio 8 = F/I of the driving force F to the 
dissipative integral 

I = 21n-^, (37) 
bh m 

where a = (3/2 V) 1 ^ 2 is the radius of a static parabolic droplet with the profile h = |a[l — (x/a) 2 ] 
and b ~ 2.082 is a constant. The driving force F = F rec — F adv is expressed through the 
equilibrium contact angles by separating the contributions of the two menisci F rec , F adv : 

F = - r \h- h m )^dx = Li^dx = F rec - F adv , (38) 

J —oo dX J —oo (XX 

where, after replacing the integration variable and extending integration to infinity in a thick mid- 
dle part of the droplet where the disjoining potential is negligible, 

precadv = / ^ec,adv(^ dh = jr^oM _ (39) 

This yields (with 9 rec = 1) 

6 = 4\n(2a/bh m y (40) 
This result is compared to the asymptotic theory in Fig|7J 

The presented asymptotic theory is based on (i) a separation into micro- and mesoscale, i.e. it 
is not valid for V too small (s becomes too small) and (ii) an expansion in S, i.e. it is not valid for 
5 too large. Assuming a precursor film of 1 nm, V = 10 6 corresponds roughly to droplets of 1 /im 
height, implying that the asymptotics is valid in the realm of microfluidics, but less so for nanoflu- 
idics. However, because for nanodroplets the micro- and mesoscale are not well separated they 
can be treated with numerical methods. It is convenient to calculate stationary moving droplets 
using continuation techniques'^ as shown, for instance, for nanodroplets moving under the influ- 
ence of a body force 49 - 50 and chemically driven droplets 41 - 42 . However, the numerical calculation 
becomes very tedious for larger drops because of the separation of scales. 

In Fig.|7]we present a comparison of asymptotic results (stretched down to V — 10 6 ) obtained 
from Eqs. (l33l) . (l34b and (l36l) . small-5 results given by Eq. d40l . and numerical continuation results 
(stretched up to V = 10 s ) for Eq. (|2TTl . For small driving 9® dv > 0.8 the overall agreement of the 
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FIG. 7: Comparison of asymptotic results (thick lines) obtained from Eqs. d33l . d34t and ( Bol l, numerical 
continuation results (corresponding thin lines) for Eq. d2Tl . and the small-5 results given by Eq. d40l (dashed 
lines). The dependence of the droplet velocity 5, on the imposed equilibrium contact angle at the advancing 
side 0ad v < rec = 1 is shown for different given droplet volumes V, as specified in the legend. 

three methods is reasonably good. For V = 10 6 the maximal deviation is below 15%, and for 
V = 10 8 it is about 5%. As expected, for larger driving 9 adv < 0.8 the results start to deviate, 
the numerical solutions of the full Eq. (OTT) give a lower velocity than the asymptotics, and more 
so for smaller 0® dv . For larger droplets this deviation starts at larger 9 adv (smaller driving). There 
are various small factors that may contribute to the deviations at small driving: (i) for V = 10 8 
the equilibrium contact angle still differs from the asymptotic value of one by about 0.3%; (ii) 
for moving droplets the precursor film thickness depends weakly on the dynamic a 49iS1 implying a 
droplet volume that is not exactly constant with changing velocity. For V = 10 8 and 9 adv = 0.5 
the precursor film thickness is about 1 .0025, i.e. for the used domain size of 10 6 the relative change 
in droplet volume is negligible (AV/V ~ 10~ 5 ). 

Surprisingly, the simple results obtained for small 5 in Ref. Q as the ratio of the driving force 
and the dissipative integral [our Eq. d40l l seem to fit the numerical data better than the asymptotic 
theory. This results, apparently, from the cancelation of different approximations. The assumed 
velocity-independent parabolic droplet shape underestimates, for instance, the dissipation at the 
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receding contact line and in the bulk, but overestimates the dissipation at the advancing contact 
line. The advantage of the asymptotic theory can be better appreciated comparing the profiles of 
the moving droplets (Fig [5]>. 

The numerical results obtained by continuation (solid lines) strongly differ from the static 
droplet shapes (dotted lines) that are the basis for the small-<5 approximation Eq. (BOl . The asymp- 
totic mesoscopic profiles [Eqs ifTTT) and (IT2b with the parameters obtained from Eqs. d33l) . d34l) and 
d36l) 1 approach the numerical results reasonably good for weak driving 9® dv = 0.8, independently 
of whether one compares profiles for identical velocity or driving (Fig|Hl(a)). For larger driving, 
the comparison of profiles for identical velocities gives better results. In general, the receding 
part is described quite perfectly. The advancing part differs because the matching is based on the 
advancing equilibrium contact angle that is smaller than the dynamical one. 

Our treatment has made it clear that the characteristics of the moving droplets depend in a 
crucial way on the kind of driving used. The droplet may be driven by body forces, as for instance, 
gravitation or Marangoni forces. In lubrication theory the latter also takes the form of a body 
force althought physically it acts at the free surface only. The driving is top-down because the 
force is fed into the system on the macroscopic scale and causes motion on all scales down to the 
microscale. One of our main results is that this type of driving cannot be described by the present 
theory because the balance of the viscous term and the capillary term in Eq. © does not account 
for the driving force. Specifically, it is not possible to use the solution of Eq. © in terms of Airy 
functions to describe droplets sliding down an incline driven by gravity. This is already obvious 
from the fact that for gravity-driven drops the advancing dynamic contact angle is larger than the 
receding on e - 49 ' 50 contrary to the characteristics of the mesoscopic solution given by Eqs. (fTTTl and 
<d. 

On the contrary, driving the droplets by a wettability gradient is bottom-up because the force 
is fed into the system on the microscale and causes motion up to the mesoscale (in our terms, no 
macroscale exists in this case because the macroscale is defined by the scale of the body forces 
that are absent by definition of the problem). 
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FIG. 8: Profiles of moving droplets for V = 10 8 as obtained by numerical continuation of stationary 
solutions [Eq. d2TT i. solid lines] are compared to the solutions of the mesoscopic asymptotic Eq. For 
the latter profiles are given that have the same 9$ (dot-dashed lines) or the same 5 (dashed line) as the 
numerical solution. Parameters are (a) numerical: 9$ = 0.8, 5 = 0.009; asymptotic (dot-dashed) 9$ = 
0.8, 5 = 0.010, s = 6.35, K/V 1 / 2 = 2.63; asymptotic (dashed) 0g d " = 0.82, 5 = 0.009, s = 2.71, 
K/V 1 / 2 = 2.71; and (b) numerical: 6f v = 0.6, 5 = 0.015; asymptotic (dot-dashed) 9f v = 0.6, 8 = 
0.019, s = 3.08, K/V 1 ' 2 = 1.98; asymptotic (dashed) 6% dv = 0.69, 5 = 0.015, s = 4.16, K/V 1 / 2 = 
2.21. For comparison the equilibrium profile of a droplet on a homogeneous substrate without wettability 
gradient is also shown (dotted lines). 
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